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Abstract. In the matrix product states approach to n species diffusion processes the 
stationary probability distribution is expressed as a matrix product state with respect to 
a quadratic algebra determined by the dynamics of the process. The quadratic algebra 
defines a noncommutative space with a SUq{n) quantum group action as its symmetry. 
Boundary processes amount to the appearance of parameter dependent linear terms in the 
algebraic relations and lead to a reduction of the SUq{n) symmetry. We argue that the 
boundary operators of the asymmetric simple exclusion process generate a tridiagonal algeb- 
ra whose irriducible representations are expressed in terms of the Askey- Wilson polynomials. 
The Askey- Wilson algebra arises as a symmetry of the boundary problem and allows to solve 
the model exactly. 
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1 Introduction 

Stochastic reaction-diffusion processes are of both theoretical and experimental interest not only 
because they describe various mechanisms in physics and chemistry [1] but they also provide 
a way of modelling phenomena like traffic flow [2], kinetics of biopolimerization [3], interface 
growth d]. 

A Markov process can be described in terms of a master equation for the probability distri- 
bution P{si, t) of a stochastic variable Si = 0, 1, 2, . . . , n — 1 at a site i = 1,2, ... ,L of a linear 
chain. A state on the lattice at a time t is determined by the occupation numbers Sj and a tran- 
sition to another configuration s • during an infinitesimal time step dt is given by the probability 
r(s, s')dt. The rates F = F*^, i, j,k,l = 0,1,2, ... ,n — I are assumed to be independent from the 
position in the bulk. At the boundaries, i.e. sites 1 and L additional processes can take place 
with rates and i?^. Due to probability conservation 

r{.s,s) = -Y,r{s',s). 

The master equation for the time evolution of a stochastic system 

^^ = ^ns,.s')Pis',t) 

s' 

is mapped to a Schrodinger equation for a quantum Hamiltonian in imaginary time 
dP{t) 



dt 



-HP{t), 



*This paper is a contribution to the Proceedings of the O'Raifeartaigh Symposium on Non-Perturbative and 
Symmetry Methods in Field Theory (June 22-24, 2006, Budapest, Hungary). The full collection is available at 
|http://www.emis.de/journals/SIGMA/LQR2006.html, 



2 



B. Aneva 



where 

j 

The probabihty distribution thus becomes a state vector in the configuration space of the quan- 
tum spin chain and the ground state of the Hamiltonian, in general non-Hermitian, corresponds 
to the steady state of the stochastic dynamics where ah probabihties are stationary. The map- 
ping provides a connection with integrable quantum spin chains and ahows for exact results of 
the stochastic dynamics with the formalism of quantum mechanics. 

Examples are the diffusion processes of particles hopping with hard core repulsion between 
lattice sites i, j with probability rates gij and subjected to reflecting (or free) boundary condi- 
tions [5]. In the symmetric case gij = gji the stochastic Hamiltonian is the SU{2) symmetric 
spin 1/2 isotropic Heisenberg ferromagnet 

H = -l + afa] + afa| - 1). 

i 

The SU{2) symmetry, yet unrevealed in the original master equation becomes manifest through 
the mapping and allows for exact results of the stochastic dynamics. The asymmetric exclusion 
process is the driven diffusive lattice gas of particles hopping with rates = <? / 1 and is 

mapped to a S'[/q(2)-symmetric XX Z chain with anisotropy A = ^^^^ — -. 



2 Matrix product state approach to diffusion models 

An algebraic approach, developed for the study of the steady state properties of Markov pro- 
cesses and related to the quantum Hamiltonian picture is the matrix-product-state ansatz [6l • 
The idea is that the stationary probability distribution, i.e. the ground state of the quantum 
Hamiltonian with nearest-neighbour interaction in the bulk and single-site boundary terms is 
expressed as a product of (or a trace over) matrices that form a representation of a quadratic 
algebra 

TfiDiDk = xiDj - XjDi, i,j,k = 0,l,...,n-l, 

determined by the dynamics of the process. 

In the following we consider diffusion processes with n species on a chain of L sites with 
nearest-neighbour interaction with exclusion, i.e. a site can be either empty or occupied by 
a particle of a given type. In the set of occupation numbers {si, S2, ■ ■ ■ , sl) specifying a configu- 
ration of the system Sj = if a site i is empty, Sj = 1 if there is a first-type particle at a site 
i, . . . , Si = n — 1 if there is an (n — l)th-type particle at a site i. On successive sites the species i 
and k exchange places with probability gikdt, where i, fc = 0, 1, 2, . . . , n — 1. With i < k, go. are 
the probability rates of hopping to the left, and gki to the right. The event of exchange happens 
if out of two adjacent sites one is a vacancy and the other is occupied by a particle, or each of 
the sites is occupied by a particle of a different type. The n-species symmetric simple exclusion 
process is known as the lattice gas model of particle hopping between nearest-neighbour sites 
with a constant rate gi^ = gki = 9- The n-species asymmetric simple exclusion process with 
hopping in a preferred direction is the driven diffusive lattice gas of particles moving under the 
action of an external field. The process is totally asymmetric if all jumps occur in one direction 
only, and partially asymmetric if there is a different non-zero probability of both left and right 

n-l 

hopping. The number of particles of each species in the bulk is conserved ^ = L and this 

1=0 
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is the case of periodic boundary conditions. In the case of open systems, the lattice gas is coupled 
to external reservoirs of particles of fixed density. In most studied examples [TJ [H] one considers 
phase transitions inducing [9] boundary processes when a particle of type k, k = 1, 2, . . . , n — 1 
is added with a rate and/or removed with a rate Lq at the left end of the chain, and it is 
removed with a rate Rq and/or added with a rate at the right end of the chain. 

For diffusion processes the transition rate matrix becomes simply Vj^- = gik and the n-species 
diffusion algebra [10] has the form 



gikDiDk - QkiDkDi = XkDi - XiD^, (2.1) 
where gik and gki are positive (or zero) probability rates, i,k = 0,1, ... ,n — 1 and Xi satisfy 



n-l 



J^x, = 0. 



i=0 

(No summation over repeated indices in equation ()2.ip .) 

The quadratic algebra has a Fock representation in an auxiliary Hilbert space where the 
n generators act as operators. For systems with periodic boundary conditions, the stationary 
probability distribution is related to the expression 

P(si, . . . , Si) = Tr {Ds,D,, ■■■Ds^). (2.2) 

When boundary processes are considered the stationary probability distribution is related to 
a matrix element in the auxiliary vector space 

P{si,...,sl) = {w\D,,D,,---Ds,\v) (2.3) 
with respect to the vectors \v) and {w\, defined by the boundary conditions 

{w\{L^Dk + X,) = 0, {R'iDk - x,)\v) = 0, 
where the x sum up to zero, because of the form of the boundary rate matrices 



L-l L-l 

R 

J- 



L\--^L], R\--Y^R 

j=0 j=0 



These relations simply mean that one associates with an occupation number Si at position i 
a matrix Dg^ = (i = 1, 2, . . . , L; A; = 0, 1, . . . , n — 1) if a site i is occupied by a /c-type particle. 
The number of all possible configurations of an n-species stochastic system on a chain of L 
sites is and this is the dimension in the configuration space of the stationary probability 
distribution as a state vector. Each component of this vector, i.e. the (unnormalized) steady- 
state weight of a given configuration, is a trace or an expectation value in the auxiliary space 
given by (j2.2p or (12. 3p . The quadratic algebra reduces the number of independent components 
to only monomials symmetrized upon using the relations (|2.ip . 

The algebra admits an involution, hence the matrices Di can be chosen Hermitian 

A = Df, gf^ = gki, Xi = -xf 

(or Di = -Df, {igik = gti)- 

The relations (|2.1|) allow an ordering of the elements D^. Monomials of given order are the 
Poincare-Birkhoff-Witt (PBW) basis in the algebra. The quadratic algebra generated by the n 
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elements obeying the n{n — l)/2 relations (|2.ip is an associative algebra with a unit e. The 
ordered monomials 

T^ni T^n2 . . . ]~)ni 

where si < S2 < • • • < s;, Z > 1 and ni,n2, ■ ■ ■ ,ni are non-negative integers, are a linear basis in 
the algebra. 

Thus, to find the stationary probability distribution one has to compute traces or matrix 
elements with respect to the vectors \v) and {w\ of monomials of the form D'^^D'^^ ■ ■ ■ D"^ . The 
problem to be solved is twofold - to find a representation of the matrices D that is a solution of 
the quadratic algebra and match the algebraic solution with the boundary conditions. Finding 
a representation of the quadratic algebra and solving the corresponding boundary problem allows 
one to conclude about the symmetries of the model. 

In the case of symmetric diffusion with all Xi distinct the diffusion algebra is of Lie-algebra 
type. The n generators Di, and e can be mapped [11] to the generators Jjk of SU{n) x U{1) 
and the mapping is invertible. The universal enveloping algebra generated by Di belongs to the 
UEA of the Lie-algebra of SU{n) x U{1). 

In the known example of exactly solvable 2- and 3-species models [TJ \T1[ [T3] of asymmetric 
diffusion through the matrix product ansatz, the solution of the quadratic algebra is provided 
by a deformed bosonic oscillator algebra, if both gik and gki differ from zero, or by infinite- 
dimensional matrices, if gik = 0. In the general n case, because of the ordering procedure, the 
solution of the quadratic algebra has to be consistent with the diamond lemma in ring theory, 
also known as the braid associativity condition in quantum groups. As shown in [10\ [TT] if 
all parameters Xi on the RHS of equation (12. ip are equal to zero, the homogeneous quadratic 
algebra defines a multiparameter quantized non-commutative space realized equivalently as a 
(^-deformed Heisenberg algebra [14^ I15j of n oscillators depending on n(n — l)/2 + 1 parameters 
(in general onn(n— l)/2 + n parameters): 

ttiaf - rittfai = 1, 
OjOj - qjiajUi = 0, 

where i < j; i,j = 0, 1, . . . , n — 1, the deformation parameters ri,qij are model-dependent 
parameters given in terms of the probability rates. The associative algebra generated by the 
elements Di in this case belongs to the universal enveloping algebra of the multiparameter 
deformed Heisenberg algebra to which a consistent multiparameter SUq{n) quantization cor- 
responds. Thus the bulk symmetry of a process without boundary conditions is SUq{n) and 
the process is equivalent to an integrable spin chain. For a non-homogeneous algebra with x- 
terms on the RHS of (12. ip . only then is braid associativity satisfied if, out of the coefficients 
Xj, Xfc, xi corresponding to an ordered triple DiD^Di, either one coefficient x is zero or two 
coefficients x are zero, and the rates are respectively related. The diffusion algebras in this 
case can be obtained by either a change of basis in the n-dimensional non-commutative space 
or by a suitable change of basis of the lower-dimensional quantum space realized equivalently 
as a lower-dimensional deformed Heisenberg algebra. The appearance of the non-zero linear 
terms in the RHS of the quantum plane relations due to boundary processes leads to a lower- 
dimensional non-commutative space and a reduction of the SUq{n) invariance. 

3 The open asymmetric exclusion process 

We consider now the two-species partially asymmetric simple exclusion process with incoming 
and outgoing particles at both boundaries. We simplify the notations, namely, at the left bounda- 
ry a particle can be added with probability adt and removed with probability 'jdt, and at the 



a+aj - qjiajaj = 0, 
aio^ - Qji^a'^CLi = 0, 
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right boundary it can be removed with probabihty (3dt and added with probabihty 5dt. The 
system is described by the configuration set si, S2, ■ ■ ■ ,sl where Sj = if a site i = 1,2, ... ,L 
is empty and Sj = 1 if a site i is occupied by a particle. The particles hop with a probability 
goidt to the left and with a probability giodt to the right, where without loss of generality we 
can choose the right probability rate 510 = 1 and the left probability rate ^foi = Q- The model 
depends on five parameters - the bulk probability rate q and the four boundary rates. The 
totally asymmetric process corresponds to q = 0. The quadratic algebra 

DiDq - qDoDi = xiDq - DiXq, xo + xi = 

is solved [16] by a pair of deformed oscillators a, a"*" (and xq = —xi = 1) 

Do = + , Di = + 



l-q VI - 9' VI - Q 

The boundary conditions have the form 

{(iDi-5Do)\v) = \v), {w\{aDo-^Di) = {w\. (3.1) 

For a given configuration (si, S2, . . . , sl) the stationary probability is given by the expectation 
value 



P{s) 



{w\Ds^Ds2 ■ ■ ■ Dsl\v) 
Zl 



where Ds^ = Di if a site i = 1,2,... ,L is occupied and Dg. = Dq if a site i is empty and 
Zl = {w\{Dq + Di)^\v) is the normalization factor to the stationary probability distribution. 
Once the representation of the diffusion algebra and the boundary vectors are known, one can 
evaluate all the physical quantities such as the current J through a bond between site i and site 
z + 1, which has a very simple form 

the mean density at a site i 



Zl 

the two-point correlation function 

{w\{Do + D^y-^DiiDo + D^y-^D^iDo + Di)^'^\ 



{s^Sj 



Zl 



and higher correlation functions. 

The matrix-product algebraic relations provide solvable (and representation independent) 
recursions for the stationary probability, the current and the correlation functions. These re- 
cursions have been obtained in earlier works [SlllZjj however they were not readily generalized 
to other models. The algebraic matrix-product method allows for generalizations to the multi- 
species case [7] and was also to the dynamical process [T5] . 

To solve the boundary problem for a process with only incoming particles at the left boundary 
and only outgoing particles at the right one (5 = 7 = in (j3.ip ) we choose the vector \v) to be 
the (unnormalized!) eigenvector of the annihilation operator a for a real value of the parameter 
V and the vector {w\ to be the eigenvector (unnormalized and different from the conjugated one) 
of the creation operator for the real parameter w: 



e,^" e-^O), (u;| = (0|e-%^'^ , (3.2) 
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where = ^ and [n] = . The factor Cq ^ in (j3.2p is due to the condition {w\v) = 1, 
which is a convenient choice in physical apphcations. According to the algebraic solution, these 
are also eigenvectors of the shifted operators Dq, Di with the corresponding relations of the 
eigenvalues 

11 w 11 V 

a 1-q q' f3 1-q \/l - q' 

The choice of the boundary vectors to be deformed coherent states |19| [20] of the g-deformed 
oscillator |20[ [2T] set underlying the algebraic solution simplifies the procedure of obtaining 
all the relevant physical quantities. One has to merely normally order all the expressions - 
monomials in a, a"*" and use the coherent states properties. The representation has the proper 
q = limit which allows to obtain all the physical quantities of the totally asymmetric process 
as well. The deformed coherent state solution in the case of only incoming particle to the left 
and only outgoing particle to the right provides the most simple and convenient approach to 
a unified description of both the partially and the totally asymmetric processes. 

The asymmetric exclusion process was studied and analyzed in relation to orthogonal polyno- 
mials technique. In the case of only injected particles at the left boundary and only removed at 
the right one the exact solution was related to the g-Hermite [16] and Al-Salam-Chihara poly- 
nomials [22]. In the general case of injected and removed particles at both boundaries with four 
nonzero boundary parameters the exact solution was shown to be related to the Askey- Wilson 
polynomials [23]. Our analysis aims to further stress the relation to the orthogonal polynomials. 
We put an emphasis on the important characteristic of the open asymmetric exclusion process 
that the bulk properties in the stationary state strongly depend on the boundary parameters. 
We find that the boundary operators generate the Askey- Wilson tridiagonal algebra whose irre- 
ducible modules are the Askey-Wilson polynomials. The boundary tridiagonal algebra reveals 
algebraic properties of the asymmetric open exclusion process which are deeply related to the 
Askey-Wilson polynomials and allow for the exact solution of the model. 



4 Symmetry of the open asymmetric exclusion process 

In the general case of incoming and outgoing particles at both boundaries there are four opera- 
tors PDi, —6Dq, —'jDi, aDo and one needs a pair of deformed oscillators to form two linear 
independent boundary operators acting on the dual vector spaces. To solve this problem within 
the matrix product ansatz we consider a particular form of a g-oscillator algebra, known in 
the literature as CUq{2) (or («,«), « > 0) algebra. It is generated by three elements with the 
defining commutation relations 

[iV,^±] = ±A±, [A_,A+] = uq^'+uq-'' 

and a central element 

TV _ 1-7V 

Q = A+A_ + u^—-^ . 

1 - q 

This is an associative algebra and any invertible transformation of the generators is admissible. 
It can be presented in equivalent forms by means of the transformations 

with 
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-N 



or 

with 

We consider now the two sets of operators 

ai = + at = + 

1 - (/ 
and 

1 — g ^ 

obeying 

oio^ — qafai = 1, Osa^ — q'~^a^a2 = 1. 
This representation ahows one to write the boundary operators in the equivalent form 
m - 6Do 

-XlP N/2 A , r -1 4 N/2 , f -^iP , x "A XoSq'^ 

= q ' A++Xodq A^q' +[- + Xo6q ]q zT, 

1-q \{l-qy J 1-q 1-q ^ 

aDo - -fDi (4.1) 



-N/2 . ^ilQ . , f ^179 \ -AT , a;oa , Xi^/q 



We separate the shift parts from the boundary operators. Denoting the corresponding rest 
operator parts by A and Al we write the left and right boundary operators in the form 

PDi - 5Dq = A - -, uDq - jDi = A + h T. 

I — q 1 — q~^ 1 — 1 — q~ 

Proposition 1. The operators A and A* defined by 

A = iSDi - 5Dq + h r, A = aDo - 7D1 - 

1 — q 1 — q ^ 1 — q 1 — q ^ 

and their q-commutator 

[A, A*]q = q^''^AA* - q~^''^A*A 
form a closed linear algebra 

[[A, A%, A], = -pA* -ujA-tj, [A*, [A, A%], = -p*A - u:A* - 7?*, (4.2) 
where the structure constants are given by 

-p = P6q~'xl{q'/' + q~'/r, -p* = a^Q~'xl{q'/' + q~'/'f, (4.3) 

-cv = (-XoXi)g"i ((/? + 6{q'/' - q-'/^) {c^il'^' " Q^'^'f + t) 

+ (-XoXi)g-i(gi/2 - q~'/^fQK, (4.4) 
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q-q- 



— = xiPxodq (^xoa + -^Y^^^l^ ^ + ^ + ^^^^ ^) ' 



q-q 
where 

K = {-XoXi)q-' (^a6 - ^^1/2^^-1/2)2 ) • (4-5) 

The proposition is straightforward to verify by using the representation for A and A* on the 
RHS of formula (14. ip . Relations (j4.2p are the well known Askey- Wilson relations 

A^A* - (g + q-^)AA*A + A*A^ = pA* + ujA + r], 
A*^A -{q + q-^)A*AA* + AA*^ = p*A + uA* + if 

for the shifted boundary operators A, A* . The algebra ()4.2p was first considered in the works of 
Zhedanov [24^ [25] and recently discussed in a more general framework of a tridiagonal algebra [261 
127] . This is an associative algebra with a unit generated by a (tridiagonal) pair of operators A, 
A* and defining relations 

[A, A^A* - p'AA*A + A*A'^ - j'{AA* + A* A) - pA*] = 0, 

[A*,A*^A - P'A*AA* + AA*^ - j'*{AA* + A* A) - p*A] = 0. (4.6) 

In the general case a tridiagonal pair is determined by the sequence of scalars /3', 7', 7'*, p, p* 
from a field K. Note that we keep the conventional notations P, 7, p for the scalars defining 
a tridiagonal pair in equation (j4.6p . however we use them with the index ' to avoid confusion 
with the boundary rates. Tridiagonal pairs have been classified according to the dependence on 
the scalars [26]. Examples are the q-Serre relations with (3' = q + q~^ and 7' = 7* = p = p* = 



[A, A^A* -{q + q~^)AA*A + A*A^] = 0, [A*, A*^A - {q + q~^)A*AA* + AA*^] = 
and the Dolan-Grady relations [28] with /?' = 2, 7' = 7'* = 0, p = , p* = k*^ 

[A, [A, [A, A*]]] =k^[A,A*], [A*,[A*,[A\A]]] = k*'^[A*,A]. 
Tridiagonal pairs are determined up to an affine transformation 

A^tA + c, A*^t*A* + c*, 

where t, t*, c, c* are some scalars. The affine transformation can be used to bring a tridiagonal 
pair in a canonical form with 7' = 7'* = 0. 

The boundary operators of the asymmetric exclusion process obeying the Askey-Wilson 
algebra ()4.2p form a tridiagonal pair with P' = q + q~^, 7' = 7'* = 0, and p, p* , uj = to* , 7], r]* 
given by equations (j4.3p - (j4.5p . The Askey-Wilson algebra possesses some important properties 
that allow to obtain its ladder representations, spectra, overlap functions. This fact immediately 
provides a solution to the boundary problem. We briefly touch the spectral problem following 
the argumentation of Zhedanov et al. (for details see |24[ 125] ). Let fr be an eigenvector of A 
with eigenvalue A,. 

Afj. — ^rfr- 

Then we can construct a new eigenstate 



/, = {Ag{A) + A*h{A) + Aok{A))fr, 
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where denotes the q-commutator ^*]q, and 

Afs = Xsfs. 

It follows from the algebra that fs will also be an eigenvector of ^, if for the new eigenvalue the 
quadratic relation holds 

Xl + X^,-{q + q~^)XrXs- p = 0. 

This yields for each state fr two neighbouring states whose eigenvalues are the roots of the above 
quadratic equation. In this parametrization the operator A is diagonal and the operator A* is 
tridiagonal 

Afr = Or+l/r+l + Kfr + Cr-lfr-1- 

The expressions for the spectrum and the matrix coefficient can be obtained explicitly, e.g. the 
spectrum reads (with xf = 1) 

The right boundary vector \v) can be thus related to a suitably chosen eigenstate of the shifted 
operator A. The algebra possesses a duality property. Due to the duality property the dual 
basis exists in which the operator A* is diagonal and the operator A is tridiagonal. The left 
boundary vector {w\ will be related to an eigenstate of the shifted operator A* in the dual basis. 
We have 

A*f; = A*/;, Af: = + b*f: + c:_j:_,. 

The overlap function of the two basis {s\r) = {fs\fr) is expressed in terms of the Askey-Wilson 
polynomials. To make the relation more transparent we rescale the generators 



The tridiagonal relations for the transformed operators read 

[A, A^A* -{q + q-^)AA*A + A* A^ + abcdq-\q - q-^fA*] = 0, 
[A*,A*^A -{q + q'^)A*AA* + AA*^ + {q - q-^)^A] = 0, 

where abed = Let p„ = pn{x;a,b,c,d) denote the nth Askey-Wilson polynomial 

depending on four parameters a, b, c, d 

with po = 1, X = y + and < g < 1. On the RHS of equation (|4.7p we use the standard 
notations for the basic hyper geometric series given in [30]. Then it can be proved (see [2^ for 
details) that pn is an eigenstate of the operator A 

Apn = (g" + q-^abcdq-'')pn, n = 0, 1, 2, . . . . 

The operator A*pn = xpn is tridiagonal and its matrix elements in the basis 

{po{x),pi{x),p2{x),...) (4.8) 



10 



B. Aneva 



are obtained from 

XPn = KPn+l + anPn + CnPn~l, P~l = 

which is the three term recurrence relation [22 of the Askey-Wilson polynomials. Let A denote 
the tridiagonal matrix 



A 



bo ai C2 
bi a2 ■ 



V 



7 



representing the operator A* in the basis (j4.8p . The spectrum of the boundary operators is 
obtained by the correspondingly rescaled shifts. We will summarize the results. In the repre- 
sentation vr with basis 

{po{x),pi{x),p2{x), ...f 

the operator Di — ^Dq is diagonal with diagonal eigenvalues 



An 



+ dg'^ + 1 + bd) 



(4.9) 



and the operator Dq — ^Di is tridiagonal 



vr ( - -Di 

a 



1 



1 



{bA^ + 1 + ac) 



(4.10) 



In the dual representation vr* with respect to the basis ()4.8p the operator — -^Di is diagonal 
with eigenvalues 



a: 



1 



l-q 



[aq-'^ + eg" + 1 + ac) 



(4.11) 



and the operator Di — 4 Do is tridiagonal 



^{D^ - -Do) 



{aA +l + bd) 



(4.12) 



The Askey-Wilson parameters a, b, c, d on the RHS of equations (|T9]) . (|4.10p and (|4.1ip . ()4.12p 
depend on the four boundary rates q, /3, 7, 6. 

The left and right boundary vectors can be chosen (up to a normalization) as 

{w\ = (po,0,0,...), It;) = (po,0,0,...)*. 

Then, the solutions to the eigenvalue equations determine uniquely, in the considered represen- 
tation of the boundary operators, the relation of the boundary rates to the four Askey-Wilson 
parameters. The explicit expressions read 



a = Kj 



d = 



where 



K± 



-{P + 6-{l-q))±y^{P + 6-{l-q)r-AP6 
2/3 
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, _ -(g + 7 - (1 - g)) ± V(a + 7 - (1 - qW 

K/_\_ — . 

We note the different relations compared to the corresponding ones in [23]. The boundary 
symmetry has rich algebraic properties and depending on its representation one has various 
identifications of the boundary rates dependence with the parameters of the Askey-Wilson 
polynomials. 

We thus conclude: The quadratic algebra of the bulk diffusion process define a quantum 
plane with the SUq{2) action as its symmetry. The boundary processes amount to the presence 
of linear terms in the quadratic algebra and lead to a reduction of the bulk symmetry. The 
Askey-Wilson algebra arises as a symmetry of the boundary problem. The boundary vectors 
are expressed in terms of the Askey-Wilson polynomials whose four parameters are related to 
the four boundary probability rates. 
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